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Abstract 

In cognitive radio (CR) networks, there are scenarios where the secondary (lower priority) users intend to 
communicate with each other by opportunistically utilizing the transmit spectrum originally allocated to the 
existing primary (higher priority) users. For such a scenario, a secondary user usually has to trade off between 
two conflicting goals at the same time: one is to maximize its own transmit throughput; and the other is to 
minimize the amount of interference it produces at each primary receiver. 

In this paper, we study this fundamental tradeoff from an information-theoretic perspective by characterizing 
the secondary user's channel capacity under both its own transmit-power constraint as well as a set of interference- 
power constraints each imposed at one of the primary receivers. In particular, this paper exploits multi-antennas at 
the secondary transmitter to effectively balance between spatial multiplexing for the secondary transmission and 
interference avoidance at the primary receivers. Convex optimization techniques are used to design algorithms 
for the optimal secondary transmit spatial spectrum that achieves the capacity of the secondary transmission. 
Suboptimal solutions for ease of implementation are also presented and their performances are compared with 
the optimal solution. Furthermore, algorithms developed for the single-channel transmission are also extended 
to the case of multi-channel transmission whereby the secondary user is able to achieve opportunistic spectrum 
sharing via transmit adaptations not only in space, but in time and frequency domains as well. Simulation results 
show that even under stringent interference-power constraints, substantial capacity gains are achievable for the 
secondary transmission by employing multi-antennas at the secondary transmitter. This is true even when the 
number of primary receivers exceeds that of secondary transmit antennas in a CR network, where an interesting 
"interference diversity" effect can be exploited. 
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I. Introduction 

Fixed spectrum allocation is the major spectrum allocation methodology for traditional wireless communication 
services. In particular, in order to avoid interference, different wireless services are allocated with different licensed 
bands. With the popularity of various wireless technologies, fixed spectrum allocation strategy has resulted in 
scarcity in radio spectrum, due to the fact that most of the available spectrum has been allocated. Because of 
the spectrum scarcity, one immediate consequence is that, while there are a lot of research activities dealing 
with technical issues related to the fourth-generation (4G) cellular systems, it is still unclear which frequency 
band is available for such systems. On the other hand, recent measurements by FCC and others have shown 
that more than 70% of the allocated spectrum in United States is not utilized [1]. Furthermore, the spectrum 
utihzation varies in space, time and frequency. This motivates the invention of cognitive radio (CR) network [2]- 
[3], which supports opportunistic spectrum sharing by allowing the secondary (lower priority) users to share the 
radio spectrum originally allocated to the primary (higher priority) users. By doing so, the utilization efficiency 
of the radio spectrum can be significantly enhanced. 

In CR networks, the primary users have a higher priority than the secondary users in utilizing spectrum resource, 
therefore, one fundamental challenge by introducing opportunistic spectrum sharing is to ensure the quahty-of- 
service (QoS) of the primary users while maximizing the achievable throughput of the secondary users. Recently, 
there has been a great deal of research related to this interesting problem. When the primary users are legacy 
systems that do not actively participate in transmit power control, the QoS of the primary users is maintained by 
introducing interference-power constraint measured at the primary receiver. That is to say, the interference power 
received at the primary receiver should be less than a threshold. Under this setup, a seminal work in [4] studied 
the channel capacity of a single secondary transmission under a set of receiver-side power constraints instead 
of the conventional transmitter-side power constraints. In [5], the ergodic capacity of the secondary transmission 
link in a fading environment is studied under instantaneous or average interference-power constraint at a single 
primary receiver. When multiple secondary users share a single-frequency channel advertised by the primary user, 
the authors in [6]- [7] utiUze the game theory to maximize the sum of the utihty functions of the secondary users 
under the interference-power constraint at some measured point, while in [8] -[9] transmit resource allocation for 
the secondary users is studied by applying the graph-theoretic models. Recent information-theoretic studies on 
CR networks can also be found in, e.g., [10]- [15]. 



Most prior research on radio resource allocation for CR networks focuses on time and/or frequency domains, 
while assuming single anteima employed at both primary and secondary transceivers. Wireless transmissions via 
multiple transmit antennas and multiple receive antennas, or the so-called multiple-input multiple-output (MIMO) 
transmissions, have received considerable attention during the past decade. Multi-antennas can be utilized to 
achieve many desirable functions for wireless transmissions, such as folded capacity increase without bandwidth 
expansion (e.g., [16]-[18]), dramatic enhancement of transmission reliability via space-time coding (e.g., [19]- 
[20]), and effective co-channel interference suppression for multiuser transmissions (e.g., [21]). However, the 
role of multi-anteimas in a CR network is yet completely understood. Generally speaking, multi-anteimas can 
be used to allocate transmit dimensions in space and hence provide the secondary transmitter in a CR network 
more degrees of freedom in space in addition to time and frequency so as to balance between maximizing its 
own transmit rate and minimizing the interference powers at the primary receivers. This motivates the research 
of this paper to be done, with an aim to address fundamentally the MIMO channel capacity of a secondary user 
under optimum spectrum sharing in a CR network. 

The main contributions of this paper are summarized as follows: 

• This paper formulates the design of capacity-achieving transmit spatial spectrum for a single secondary link 
in a CR network under both its own transmit-power constraint and a set of interference-power constraints 
at the primary receivers as a sequence of convex optimization problems. Thanks to convexity of these 
formulated problems, efficient numerical algorithms are proposed to obtain the optimal secondary transmit 
spatial spectrum for any arbitrary number of secondary transmit and receive anteimas, as well as primary 
receivers each having single or multiple antennas. 

• In the case where the secondary user's channel is multiple-input single-output (MISO), i.e., there is only a 
single anteima at the secondary receiver, this paper proves that beamforming is the optimal strategy for the 
secondary transmission. For the special case where there is only one single-antenna primary receiver, we 
are able to derive the closed-form solution for the optimal beamforming vector at the secondary transmitter. 
For the more general MIMO case where multiple antennas are equipped at both the secondary transmitter 
and receiver, this paper presents two suboptimal algorithms to tradeoff between spatial multiplexing for the 
secondary transmission and interference avoidance at the primary receivers. One algorithm is based on the 
singular- value decomposition (SVD) of the secondary user's MIMO chaimel directly and is thus referred 
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to as the Direct-Channel SVD (D-SVD); and the other is based on the SVD of the secondary MIMO 
channel after the projection into the null space of the channel from the secondary transmitter to the primary 
receivers (thereby removing completely the interference at all primary receivers) and is thus referred to as 
the Projected-Channel SVD (P-SVD). 

• In the case of multiple primary receivers with single or multiple antennas, a hybrid D-SVD/P-SVD algorithm 
is proposed to remove partially, as opposed to completely in the case of P-SVD, the interferences at the 
primary receivers. An interesting and novel interference diversity effect is discovered and exploited. 

• At last, this paper extends the developed algorithms for the single-channel transmission to the case of multi- 
channel transmission whereby the secondary transmitter is able to adapt transmit resources like transmit 
spatial spectrum, power and rate in both space and tune/freqeuncy for opportunistic spectrum sharing. This 
paper shows that the multi-channel resource allocation problem can be efficiently solved using the Lagrange 
dual-decomposition method that decomposes the original problem into a set of smaller-size subproblems 
each for one of the sub-channels. 

In [4], Gastpar also considered the MIMO channel capacity for the spectrum sharing scenario. Since only the 
interference-power constraint is considered in [4], it requires that there exist multiple primary receivers/antennas 
such that the channel matrix from the secondary transmitter to the primary receivers/antennas is invertible, and 
by doing so, an upper-bound of the secondary user's MIMO channel capacity is derived by transforming the set 
of receiver-side power constraints into a transmitter-side power constraint. In contrast, in this paper we consider 
the MIMO channel capacity for the secondary user under both the interference-power constraints at the primary 
receivers as well as an explicit transmit-power constraint for the secondary user. With addition of this explicit 
transmitter-side power constraint, we are able to quantify the exact channel capacity for the secondary user for 
any arbitrary number of primary receivers/antennas. 

This paper is organized as follows. Section |II] provides the system model of a CR network under opportunistic 
spectrum sharing, and presents the general problem formulation for transmit optimization of the secondary user. 
Section |lll] presents the solution for the simplest case where there is only one single-antenna primary receiver, 
and both primary and secondary users share the same single-channel for transmission. Section |IV] and Section 
rvl extend the solutions to incorporate multiple primary receivers/antennas and the multi-channel transmission, 
respectively. Section IVTl provides the simulation results. Finally, Section IVlIl concludes the paper. 
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The following notations are used in this paper. \S\ denotes the determinant, and Tr(S') the trace of a square 
matrix S, and S ^ means that S is a positive semi-definite matrix. For any general matrix M, denotes 
its conjugate transpose, and Rank(Af) denotes its rank. I denotes the identity matrix. E[-] denotes statistical 
expectation. C^^v denotes the space of x x y matrices with complex entries. The distribution of a circularly- 
symmetric-complex-Gaussian (CSCG) vector with the mean vector x and the covariance matrix S is denoted 
by CJ\f{x,'S), and ~ means "distributed as". The quantity m.m{x,y) and max{x,y) denote, respectively, the 
minimum and the maximum between two real numbers, x and y, and = max(x,0). The quantity ||a;|| 

denotes the Euclidean norm of a vector (also for a scalar) x, i.e., ||a;|p = Tt^xx"^). Re(x) and Im(x) denote the 
real and imaginary part of a complex number x, respectively. 

II. System Model and Problem Formulation 

This paper considers a CR network with K primary receivers and a single pair of secondary transmitter and 
receiver as shown in Fig. [T] It is assumed that all the primary users and the secondary user share the same 
bandwidth for transmission. This paper considers the scenario where multiple antennas are equipped at the 
secondary transmitter and possibly at the secondary receiver and each of the primary receivers. It is assumed that 
the MIMO/MISO channels from the secondary transmitter to the secondary and primary receivers are perfectly 
known at the secondary transmitter. Under such assumptions, the secondary transmitter is able to adapt its transmit 
resources such as transmit rate, power, and spatial spectrum based upon the channel knowledge so as to optimally 
balance between maximizing its own transmit throughput and avoiding interferences at the primary receivers. In 
practice, the channel from the secondary transmitter to the primary receiver can be obtained at the secondary 
transmitter by, e.g., periodically sensing the transmitted signal from the primary receiver provided that time- 
division-duplex (TDD) is employed by the primary transmission. In a fading environment, there are cases where 
it is difficult for the secondary transmitter to perfectly estimate instantaneous channels. In such cases, the results 
obtained in this paper provide capacity upper-bounds for the secondary transmission in a CR network. 

Consider first a single-channel transmission (e.g., narrow-band transmission with deterministic channels) for 
both primary and secondary users. The extension to the multi-channel transmission (e.g., multi-tone transmission 
in frequency or consecutive block-fading channels in time) is considered later in Section |V] of this paper. In the 
single-channel case, the secondary transmission can be represented by 

y(n) = Hx{n) + z{n). (1) 



In the above, H G £^Mr:sxMt,s denotes the secondary user's channel (assumed to be full rank) where M^^s and 
s are the number of antennas at the secondary receiver and transmitter, respectively; y{n) and x{n) are the 



index; z{n) is the additive noise vector 
Let the transmit covariance matrix (or 



received and transmitted signal vector, respectively, and n is the symbo 
at the secondary receiver, and it is assumed that z{n) ~ CAA(0, 
spatial spectrum) of the secondary user be denoted hy S, S = ¥,[x{n)x^ {n)], where the expectation is taken 
over the code-book. Since this paper characterizes the information-theoretic limit of the secondary transmission, 
it is assumed that the ideal Gaussian code-book with infinitely large number of codeword symbols is used, i.e., 
x{n) CM{0, S),n = 1,2, ... ,00. The transmit covariance matrix S can be further represented by its eigenvalue 
decomposition expressed as 

S = VSV^ (2) 

where V G C*^* ^^'^, V'^V = I, contains the eigenvectors of S, and is also termed in practice as the precoding 
matrix because each column of V is the precoding vector for one transmitted data stream; d, d < Mt^s, is usually 
referred to as the degree of spatial multiplexing because it measures the number of transmit dimensions (or 
equivalently, the number of data streams) in the spatial domain; S is a d x d diagonal matrix, and its diagonal 
elements, denoted by cji, 02, ■ ■ ■ , (Jd, are the positive eigenvalues of S, and also represent the assigned transmit 
powers for their corresponding data streams. Notice that Rank(S') = d. If d = I, the corresponding transmission 
strategy is usually termed as beamforming while in the case of d > 1, it is termed as spatial multiplexing. The 
transmit power P for each block is limited by the secondary user's own transmit power constraint denoted by 
Pt, i.e., it holds that P = Tr(5) = Ya=i ^"i < Pt- 

Assuming that there are K primary receivers in the CR network, each equipped with receive antennas, 
k = 1,. . . ,K. For each primary receiver, there may be a total interference-power constraint over all receive 
antennas or a set of interference-power constraints applied to each individual receive antenna. The former case 
can be expressed as 

Y,9k,jSgl,<'^k, k = l,...,K, (3) 

where g/^j € C^^^^' = represents the channel from the secondary transmitter to the j-th receive antenna of the A;-th 

primary receiver, and Ffc is the total interference-power constraint over all receive antennas for the A;-th primary 

'The noise at the secondary receiver may also contain the interference from the primary transmitters (not shown in Fig. [TJ and is thus 
non-white in general. However, by applying a noise-whitening filter at the secondary receiver and incorporating this filter matrix into the 
channel matrix H, the equivalent noise at the secondary receiver can be assumed to be approximately white Gaussian. 
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receiver. Let Gk G ([^MkxMt.s (assumed to be full-rank) be the equivalent channel matrix from the secondary 
transmitter to the k-th primary receiver obtained by stacking all gi^j, j = 1, . . . ,Mk, into Gk- Using G^'s, (|3]) 
can be thus rewritten as 

Tr(GkSGl^ <rk, k = l,...,K. (4) 

The latter case can be represented as 

9k,jSg\j<lk, j = l,...,Mk,k=l,...,K, (5) 

where 7^ is the interference-power constraint applied to each antenna of the k-th primary receiver, and is assumed 
to be identical for all of its receive antennas. Notice that if = M^^k, the per-antenna power constraint in ^ 
is more stringent than the total-power constraint in ©. On the other hand, dSjl can be considered as a special 
case of (O because if each receive antenna is treated as an independent primary receiver, dSjl is equivalent to (|3]l 
if a total number of M^^p = XIa^i single-antenna primary receivers are assumed. Therefore, for the rest of 
this paper, the total interference-power constraint in Q or (|4]l is considered as the generalized interference-power 
constraint at each primary receiver. 

The use of total interference-power constraint at each primary receiver can be further justified by the following 
theorem: 

Theorem 1: The capacity loss of the A:-th primary transmission. A: = 1, 2, . . . , K, due to the interference from 
the secondary transmitter with transmit covariance matrix S that satisfies the set of total interference-power 
constraints in (|4]| is upper-bounded by min(Mfc, A'^fc) log2 ^1 + bits/complex dimension, where and 
are the number of transmit and receiver antennas for the fc-th primary user, respectively, and (pk is the additive 
white Gaussian noise power at its receiver. 

Proof: Please refer to Appendix IH ■ 
From Theorem[Tl it follows that by imposing the total interference-power constraint, it is ensured that the capacity 
loss of each primary user due to the secondary transmission is well regulated. Notice that the obtained upper- 
bound for the capacity loss is not a function of the primary/secondary user's channels, or their transmit covariance 
matrices. By choosing to be sufficiently small compared with the noise power at each primary receiver, e.g., 
Ta; ^ 4>k^ the capacity loss due to the secondary transmission can be made arbitrarily small. 

At last, we formulate the main problem to be addressed in this paper. We are interested in the design of the 
spatial spectrum S at the secondary transmitter so as to maximize its transmit rate under both its own transmit- 



power constraint and a set of total interference-power constraints at K primary receivers. Accordingly, the optimal 
S can be obtained by solving the following problem (PI): 

/ + HSH^ 



Maximize log2 



(6) 



Subject to Tr(5) < Pt, (7) 
Tr(Gfc5Gl) <Tk, k = l,...,K (8) 
S^^. (9) 

The above problem maximizes the secondary user's channel capacity (in bits/complex dimension) obtained by 
computing the mutual information [22] between the channel input and output in ([T]), assuming that the secondary 
user's channel is also known perfectly at the secondary receiver. The constraints in ([7]) and ([U correspond 
to the secondary transmitter power constraint and the interference-power constraints at the primary receivers, 
respectively. ^ is due to the fact that the spectrum matrix S must be positive semi-definite. The problem at 
hand is a convex optimization problem [23] because its objective function is a concave function of S and all of 
its constraints specify a convex set of S. Therefore, for any arbitrary number of secondary transmit and receive 
antennas, and primary receivers each having single or multiple antennas, this problem can be efficiently solved 
by using standard convex optimization techniques, e.g., the interior-point method [23], for which the details 
are omitted here for brevity. In the sequel, we will investigate further into this problem so as to provide more 
insightful solutions for it that may not be obtainable from solely a numerical optimization perspective. We will 
first consider the simplest case where there is only one single-antenna primary receiver at present in Section [nil 
and then consider the general case of multiple primary receivers/antennas in Section |IVl 

III. One Single-Antenna Primary Receiver 

The scenario where there is only one single-antenna primary receiver is considered in this section. Consequently, 
the channel from the secondary transmitter to the primary receiver is MISO, and can be represented by a vector 
g G C^^^^' ". Let 7 denote the maximum interference power tolerable at the primary receiver. Problem PI can 



be then simplified as (P2) 

Maximize logs I + HSH^ (10) 

Subject to Tr(S') < Pt, (11) 

gSg^ < 7, (12) 

5^0. (13) 

A. MISO Secondary User's Channel 

Consider first the case where there is only a single antenna at the secondary receiver, i.e., M^^s = 1 and the 
secondary user's channel is also MISO. Hence, H is in fact a vector and for convenience is denoted hy H = h, 
h G C^^*^"'. In this case, we are able to derive the closed-form solution for the optimal S. First, the following 
two lemmas are needed: 

Lemma 1: In the case of MISO secondary user's channel, the optimal S for Problem P2 satisfies Rank(S) = 1. 

Proof: Please refer to Appendix |II1 ■ 
Lemma [T] indicates that in the MISO case, beamforming is indeed optimal for the secondary transmitter. 
Therefore, S can be written in the form of 5 = vv'^jV G C^^* =^^. Problem P2 can be then simplified as 

Maximize logg (l + | P) (14) 
Subject to ll^'ll^ <-Pt, (15) 
\\9v\\^<j. (16) 

The second lemma then provides the optimal structure of v: 

Lemma 2: In the case of MISO secondary user's channel, the optimal beamforming vector v is in the form 
of a^g + j3vh±, where g = and h± = j^fy^' h± = — {g'^h^)g, and /3„ are complex weights. 

Proof: Please refer to Appendix JII] ■ 

Lemma |2] states that the optimal beamforming vector should lie in the space spanned jointly by g'^ and the 
projection of into the null space of g^ . Using Lemma |2] and let h.^ = aug + /3hh±, the optimal weights ay 
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and (]y can be then obtained by considering the following equivalent problem of that in (fT4l)-(fT6l): 

Maximize W^la^ + PlfS^W'^ (17) 

Subject to \\ay\\'^ + \\Py\\'^ < Pt, (18) 

llffll'lla.ll' <7- (19) 

The above problem can be solved by using standard geometry. Notice that the solution for the above problem is 
trivial in the case of ah = or Ph = 0. Hence, without loss of generality, it is assumed that 7^ and (3h 7^ 0. 
To summarize, the following theorem is established: 

Theorem 2: In the case of MISO secondary user's channel, the optimal transmit covariance matrix S for 
Problem P2 can be written in the form of S = vv\ where v = UvQ + Pvh±, and ay and (3y are given by 

. Case I: If 7 > ^CW^p 



. Case II- If -Y < R 
, L,ase 11. II 7 < ||„^||2+p^||2^t, 



ah o fZ T~ Ph 



\g\\\Wh\y V iiffip ii/3hir 

In Theorem [2j Case I corresponds to the interference-power constraint in ([T6l ) being inactive. In this case, it 
can be verified that v = Therefore, the optimal beamforming vector is identical to that obtained by the 

pre-maximal-ratio-combining (MRC) principle for the conventional MISO point-to-point transmission. On the 
other hand. Case II corresponds to the interference-power constraint being active and, hence, the transmit power 
1 1 a„ I P allocated in the direction of g needs to be regulated by 7. 

B. MIMO Secondary User's Channel 

In the case that the secondary user's channel is MIMO, i.e., Mt^s > 1 and M^^s > 1, there is in general 
no closed-form solution for the optimal S, and Problem P2 needs to be solved numerically. Unlike the MISO 
case, it is possible that Rank(S^) > 1 in the MIMO case, which implies that spatial multiplexing is optimal 
instead of beamforming. Without the interference-power constraint in ([T2l) . the optimal S for Problem P2 can be 
obtained from the SVD of H, along with the water-filling (WF) -based power allocation (e.g., [16], [22]). The 
SVD-based transmission is not only capacity-achieving, but of more practical significance, it diagonalizes the 
secondary user's MIMO channel matrix and decomposes it into parallel AWGN sub-channels, and thereby, reduces 
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substantially the overall encoding/decoding complexity since independent encoding/decoding can be applied over 
these sub-channels. Unfortunately, when the interference-power constraint in ([T2l ) is applied, the optimal S is in 
general different from that obtained from the channel SVD, and hence does not diagonalize the secondary user's 
MIMO channel. As a result, sophisticated encoding and decoding methods need to be used in order to achieve 
the secondary user's channel capacity. One capacity-achieving method is to use a single constant-rate Gaussian 
code-book [24] to encode the information bits and then spread the coded symbols to all transmit antennas. At 
the receiver, the maximum-likelihood (ML) -based detection and iterative decoding (e.g., [25]) are applied to 
decode the whole codeword. Alternatively, the capacity can also be achieved by using multiple variable-rate 
Gaussian code-books each for one of d data streams. At the receiver, these data streams are decoded by optimum 
decision-feedback-based successive decoding (e.g., [26]-[27]). 

This paper presents two suboptimal algorithms that obtain closed-form solutions for S in the case of MIMO 
secondary user's channel. The computational complexity for both algorithms is much lower than that of the 
interior-point method for Problem P2. Furthermore, both algorithms are based on the SVD of the secondary user's 
channel matrix and decompose the secondary MIMO channel into parallel sub-channels. The main difference 
between these two algorithms lies in that, in the first one, the channel decomposition is based on the SVD of 
H directly and is thus referred to as the Direct-Channel SVD (D-SVD), while in the second one, the channel 
decomposition is based on the SVD of the projection of H into the null space of g^, and is thus referred to 
as the Projected-Channel SVD (P-SVD). Notice that these two algorithms can also be applied in the previous 
case of MISO secondary user's channel to obtain closed-form (but in general suboptimal) solutions for transmit 
beamforming vector v. 

1) Direct-Channel SVD (D-SVD): The eigenvalue decomposition of S is represented by the precoding matrix 
V and the power allocation S in Q. In the D-SVD, the precoding matrix V is obtained from the SVD of 
H, which can be expressed as H = QA^/^f/^ where Q G ^^U.^xM. ^nd U G c^^h.^xM, ^re matrices with 
othornormal columns, Mg = min(Mj .j, M^^s), and A is a Ms x Mg diagonal and positive matrix for which its 
diagonal elements are denoted by Ai > A2 > • • • > Aj\/,. Let V = U, and furthermore, y{n) = Q^y{n) and 
x{n) = Ux{n), the secondary user's MIMO channel in (dJ can be equivalently written as 

yin) = Q^HUx{n) +z{n), (20) 
= A^/2^(n) + z(n), (21) 
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where z{n) = Q'^ z{n) and z{n) ~ C7\A(0,/). Notice that the D-SVD diagonalizes the MIMO channel and 
decomposes it into Mg parallel sub-channels with channel gains \/Ai, i = 1, . . . , Ms- 

Let U = [ui,U2, . . . , wjv/J, and ai = \ \gui\\'^ ,i = 1, . . . , Ms- Considering the equivalent channel in (|2T]) and 
noticing that W\x{n)x\n)\ = XI, the optimal power assignments cjj's can be obtained by solving the following 
equivalent problem (P3) derived from Problem P2: 



Maximize log2 (1 + AjCJj) (22) 

■t=i 

Subject to ^o-j < Pt, (23) 

i=l 
A/, 

X]«i^i^7, (24) 

j=l 

(Ji>0,Vi. (25) 
The optimal cji's for the above problem can be shown to have the following form of multi-level WF solutions: 



ai= — , i = l,...,M„ (26) 

where v and /i can be shown to be the non-negative Lagrange multipliers [23] associated with the transmit-power 
constraint in (1231 ) and the interference-power constraint in (l24l ). respectively, and can be obtained by the following 
algorithm (Al): 

• Given ^ G [0, /i] 

• Initialize /J,min = 0, ^max = A 

• Repeat 

1. Set /i < 2(/^iiiin ~^ /^max)- 

2. Find the minimum u, v >0, with which — j-)^ ^ -Ft- Substitute the obtained u into (l26l) 
to obtain cjj's. 

3. Update /x by the bisection method [23]: If ajCTj > 7, set ;Umin ^ otherwise, 

Mmax ^ Z^- 

• Until /Xmax — /^min < ^fi whcrc 6fj, is a small positive constant that controls the algorithm accuracy. 
Remark 3.1: If the interference -power constraint in (l24l ) is inactive, it follows from the Karush-Kuhn-Tacker 

(KKT) conditions [23] that /x = 0. From (l26l ). the allocated powers can be then written as (^i = (j; — x") ' ^ — 
1, . . . , Ms, which become equal to the standard WF solutions with a constant water-level ^. Hence, the D-SVD 
is indeed optimal if 7 is sufficiently large such that the interference -power constraint is inactive. 
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Remark 3.2: If the secondary user's channel is MISO, it is not hard to verify that the D-SVD algorithm results 
in S = vv^, where v = Comparing it with the optimal S by Theorem |2l it follows that the 

D-SVD is optimal only if Case I of Theorem |2] is true, i.e., 7 is sufficiently large such that the interference-power 
constraint is inactive. 

2) Projected-Channel SVD (P-SVD): In the P-SVD, S is designed based on the so-called zero-forcing (ZF) 
criterion so as to completely avoid any interference at the primary receiver. It is noted that the ZF criterion has 
also been used in the MIMO broadcast channel (MIMO-BC) for design of downlink precoding that removes any 
co-channel interference between users (e.g., [28]-[29]). For the P-SVD, the secondary user's channel H is first 
projected into the null space of g"^ as 

H^ = h{i- gg^) . (27) 

Let the SVD of the projected channel H ±^ be expressed as H ±^ = Q^Ax(?7x)1^. Then, the precoding matrix 
V for S is taken as V = U _\_. From (l27l) . by multiplying both the left-hand-side (LHS) and right-hand-side 
(RHS) with g, it can be verified that (L/"_l)^S' = 0. Since S = U aT,{U it then follows that the resultant 
interference power at the primary receiver is zero, i.e., gSg^ = 0, and hence the interference -power constraint 
([T2I) for Problem P2 is always satisfied for any 7, 7 > 0. 

Like the D-SVD, the P-SVD also diagonalizes the secondary user's MIMO channel. To see it, let y{n) = 
{Q±)'^y{'n) and x{n) = U±x{n), and the secondary channel in ([T]) can be equivalently represented by 



y(n) = {Q^yHU^x{n) + z{n), (28) 

= {Q^)HH^ + Hgg^)U^x{n) + z{n), (29) 

= {Q^)^HxU±x{n)+z{n), (30) 

= {A^)^/^xin) + z{n), (31) 

where z{n) = {Qi_)'^z{n) and z{n) ~ CAA(0,/), (US is from and ^ is by g^U±_ = 0. Because the 
projection in ( |27l ) reduces the channel rank at most by one, Aj^ has = min(M( 5 — l,Mr^s) diagonal 



elements, denoted by A;f > . . . > ■ Hence, the secondary MIMO channel is decomposed into sub- 
channels with channel gains ^ ^ ~ 1) • • • The power allocation S for these sub-channels can be then 
obtained by considering the equivalent channel (|3TI ) with W\x{n)x\n)] = Xl, as the standard WF solutions 
(jj = ( I/' — p- j , i = 1, . . . , Mg, where u' is the constant water-level such that J2i=i = Pt- 
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Remark 3.3: If the interference-power constraint for Problem P2 is indeed 7 = 0, it is conjectured that the 
P-SVD is optimal. This conjecture is proved in Appendix |IVl 

Remark 3.4: If the secondary user's channel is MISO, it is not hard to verify that the P-SVD algorithm results 
in S = vv\ where v = V^'p^^x- Comparing it with the optimal S by Theorem |2l it follows that the P-SVD 
is in general suboptimal unless 7 = and then Case II of Theorem |2] applies. 

3) Performance Comparison: Comparing the computational complexity of the D-SVD and P-SVD, it is noticed 
that the former has a larger complexity for determining optimal power allocations due to the multi-level WF, 
while so does the latter for obtaining the precoding matrix due to the additional channel projection. In the 
following, the maximum achievable rate (or the capacity) for the secondary link by the optimal S, and that by 
the two SVD-based algorithms are compared under two extreme scenarios: Pt ^ and Pt — > cxd, both under 
the assumption that 7 is finite and 7 > 0. Notice that the former case can also be considered as asymptotically 
low signal-to-noise ratio (SNR) at the secondary receiver (if the transmit power P takes its maximum value Pt) 
while the latter case as asymptotically high SNR. 

As Pt 0, the D-SVD is the optimal solution. This can be easily verified similar like in Remark 13.11 by 
observing that the interference -power constraint in ([T2l) eventually becomes inactive as Pt —>■ 0. In contrast, the 
P-SVD may incur a non-negligible rate loss in this case. This is so because the ZF-based projection causes part 
of the secondary user's channel energy to be lost. Since as Pt — > 0, both SVD-based algorithms assign the total 
transmit power to the sub-channel with the largest channel gain (Ai for the D-SVD and A^*" for the P-SVD), it can 
be easily verified that the achievable rates of the D-SVD and P-SVD become log2(l + ^iPt) and log2(l + ^iPt), 
respectively. Notice that Ai > A^*". Because log(l + x) = x as x — > 0, the following theorem is obtained for 
the achievable rate by the optimal S, the D-SVD and the P-SVD, denoted by i?opt> Pd-svd> and i2p-svD> 
respectively: 

Theorem 3: As Pt 0, i?opt = ^D-SVD = and -Rp_svD = 

On the other hand, as Pt —>■ 00, the achievable rate by the P-SVD becomes close to the secondary channel 
capacity achievable by the optimal S. This is so because as Pt increases and eventually exceeds some certain 
threshold, any additional transmit power needs to be allocated into the projected channel in order not to 
violate the interference -power constraint. From (|3TI ). it can be easily verified that as Pt —>■ 00, the achievable 
rate of the P-SVD has a linear increase with log2 Pt by a factor of = min(Mr,s, Mt ,5 — 1), which is usually 
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termed as the pre-log factor, the degree of freedom, or the spatial multiplexing gain in the literature. The following 
lemma states that i?opt also has a spatial multiplexing gain of M^, indicating that the P-SVD is indeed optimal 
(in terms of the achievable spatial multiplexing gain) as Pt — > oo. 

Lemma 3: As Pt oo, j^^^ = m.m{Mr,s, Mt^s - 1)- 

Proof: Please refer to Appendix |Vl ■ 
In contrast, the D-SVD may suffer a significant rate loss when Pt oo. This can be observed from ((24l) . Suppose 
ai > 0,Vi, the total transmit power P is upper-bounded as P = — ISn"^"' regardless of the secondary 

user's own transmit power constraint P^. As a result, the achievable rate of the D-SVD eventually gets saturated 
as Pt —>■ oo, and the corresponding spatial multiplexing gain can be easily shown equal to zero. The above results 
are summarized in the following theorem: 

Theorem 4: As Pt ^ oo, = = min(M^,„ Mt,s - 1), and = if > 0, i = 1, . . . , M,. 

IV. Multiple Primary Receivers 

In this section, the general case of multiple primary receivers each having single or multiple receive antennas 
as shown in Fig. [J is studied. Hence, Problem PI in ©-(IQ]) is considered. 

A. MISO Secondary User's Channel 

First, we consider the special case where the secondary user's channel is MISO. Like Lemma [T] in the case 
of one single-antenna primary receiver, it can be shown that the optimal transmit covariance matrix S in the 
case of multiple primary receivers/antennas still remains as a rank-one matrix when the secondary user's channel 
is MISO, i.e., beamforming is optimal. Hence, it follows that S = vv'^ . However, a closed-form solution for 
the optimal beamforming vector v seems unlikely when the total number of antennas from all primary receivers 
Mj. p = X^^i is greater than one. Nevertheless, the optimal rank-one S for Problem PI can still be obtained using 
numerical optimization techniques like the interior-point method. Alternatively, the two SVD-based algorithms 
can also be modified to obtain suboptimal v, as will be illustrated later in Section HV-BI when the more general 
case of MIMO secondary user's channel is addressed. 

In this subsection, we first consider Problem PI by substituting S = vv'^ into the problem. Then, based on 
this reformulated problem, we present alternative numerical optimization techniques that, in most cases, can be 
more efficient in terms of computational complexity than the interior-point method. The reformulated problem 
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of Problem PI when the secondary user's channel is MISO, i.e., H = h, can be expressed as (P4) 

1 1 2 

Maximize H^i'll (32) 

Subject to \\v\\'^ < Pt, (33) 

\\Gkvf<rk, k = l,...,K. (34) 

Notice that for the above problem although the two constraints specify a convex set for v, the objective function 
is non-concave of v and hence renders Problem P4 non-convex in its direct form. However, an interesting 
observation here is that if v satisfies the two constraints (1331 ) and (l34l ). so does e^^v for any arbitrary 9, and the 
value of the objective function is maintained. Thus, without loss of generality we can assume that hv is a real 
number. Using this assumption. Problem P4 can be rewritten as 

Maximize Re{hv) (35) 

Subject to liii{hv) = 0, (36) 

< Pt, (37) 

\\Gkvf<Tk, k = l,...,K. (38) 

It can be shown that the above problem can be cast as a second-order cone programming (SOCP) [23], which 
can be solved by standard numerical optimization software. 

Since Problem P4 is indeed convex, and it can be further verified that the Slater's condition [23] that requires 
that the original (primal) problem needs to have a non-empty interior of the feasible set is satisfied for this 
problem, it follows that the strong duality holds for the problem at hand, which ensures that the duality gap 
between the original problem and its Lagrange dual problem is zero, i.e., solving the Lagrange dual problem 
is equivalent to solving the original problem. Motivated by this fact, as follows we will present an alternative 
solution for Problem P4 by considering its Lagrange dual problem. First, the Lagrangian [23] of Problem P4 can 
be expressed as 

K 

Civ, u, = v^h^hv - V (v^v -Pt)-Yl (v^GlGkV - Tu) , (39) 

k=l 

where v and fik, k = 1, . . . , K, are non-negative Lagrange multiplier (dual variable) associated with the transmit- 
power constraint in ([33]) . and each interference-power constraint in (l34l) . respectively. Let be the set specified 
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by the constraints (l33l) and (l34l) . the Lagrange dual function [23] can be then defined as 

f{^, {fJ'k}) = C{v, I/, {fik})- (40) 
The Lagrangian dual problem [23] is then defined as 

>n™^>nvi.-^('''^^^^^- ^^^^ 

u>0,fik>0,vk 



Notice that the Lagrangian in ( |39l ) can be rewritten as 

(h^h -vI-Y, ^fcG^Gfc^ v + uPt + Y, ^^kTk. (42) 

V fe=i / fc=i 

From the above expression, it follows that for solving the maximization problem so as to obtain the dual function 
in ( |40l ). the dual variables u and {/x^} must satisfy 

K 

h)h -uI-Y, l^kGlGk ^ 0, (43) 

k=l 

because otherwise, the optimal v that maximizes the Lagrangian in (|42l) will go unbounded to infinity, which 
contradicts the facts that the optimal value of the primal problem is bounded, and the duality gap is zero. Using 
(|43] ). the dual problem in (|4TI) can be expressed as (P5) 

K 

Minimize uPt + ^ fik^k (44) 
k=l 

K 

Subject to -h^h + vl + 'YfikGlGkhO, (45) 

A;=l 

z/ > 0, (46) 
/ifc >0, k = l,...,K. (47) 

The above problem can be efficiently solved using the standard semi-definite programming (SDP) software [23]. 
From the dual optimal solutions, the primal optimal solution for v can be obtained as any vector (may not be 
unique) that satisfies (^h^h — ul — Ylik=i l^kG\G]^ u = 0. 

The computational complexity of three different methods for solving Problem PI, namely, the interior-point 
method, the SOCP for P4, and the SDP for P5, in the case of MISO secondary user's channel are compared as 
follows. All three methods have the similar order of complexity in terms of the number of (real) variables in each 
corresponding problem. For the interior-point method, the unknown S has Mt,s real variables and 
complex variables, results in a total number of real variables equal to M^^. The SOCP has the unknown v that has 
2Mt^s real variables. The SDP has the unknown dual variables u and {jik} that in total contribute to K + 1 real 



18 

variables. Consequently, the SOCP has a much lower computational complexity than the interior-point method. 
So does the SDP than the SOCP when, e.g., K < Mt^s- At last, it is noted that the SOCP and SDP have also 
been applied for design of downlink precoding for MIMO-BC channels (e.g., [31]-[32]). 

B. MIMO Secondary User's Channel 

When the secondary user's channel is MIMO, in general we have to resort to the interior-point method for 
solving Problem PI. Alternatively, the two SVD-based solutions developed in Section |lll] in the case of one 
single-antenna primary receiver can also be modified to incorporate more than one primary receivers/antennas, 
as will be shown in this subsection. 

1) Direct-Channel SVD (D-SVD): For the D-SVD, the precoding matrix V remains to be U obtained from 
the SVD of H regardless of the number of primary receivers or antennas. However, the power allocation Xl 
needs to be adjusted in order to incorporate multiple interference-power constraints in ([3]). By introducing a set 
of non-negative Lagrange multipliers, /ii, . . . , ^k, each associated with one of K interference-power constraints 
in Like Problem P3 in the case of one single-antenna primary receiver, the optimal power assignments fjj's 
in the case of multiple primary receivers/antennas can also be obtained as the multi-level WF solutions: 

^^=( ^^K ^M. y) . ^ = 1,---,M„ (48) 

V^ + Z]fc=iz]j=i«i,fcj7^fc -^V 

where ai,fc,j = | {g^jUiW'^. Algorithm Al also needs to be modified so as to iteratively search for the optimal u and 
^fc's. Instead of updating a single fi by the bisection method in Algorithm Al, //^'s can be simultaneously updated 
by the ellipsoid method [30] by observing that 7^ — Yli^i Ylj=i ^i,j,k^i a sub-gradient of ni^, k = 1, . . . , K. 

2) Projected-Channel SVD (P-SVD): Let G e C*^'^"^^' ' denote the channel from the secondary transmitter 
to all primary receivers/antennas by taking each Qj^j as the (X^^/^i M^'-i +i)-th row of G, where Mq = 0. 
Alternatively, G = [GJ, . . . ,GJ^]'^. Let the SVD of G be expressed as G = QqAq'^Uq. In order to avoid 
completely the interference at all primary receivers/antennas, the P-SVD obtains the precoding matrix V from 
the SVD of H±, which is the projection of H into the null space of G"^ , i.e., 

H^=H{I-UgUI). (49) 

Notice that the above matrix projection is non-trivial only when Mt^s > Mr^p, because otherwise a null matrix 
H± is resulted. Let the SVD of H± be expressed as H± = Q ±{A.±)^/'^ (U . Like the case of a single- 
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antenna primary receiver, it can be verified that the secondary user's MIMO channel can be decomposed into 
m.m{Mr^s, Mt^s — ^r.p) sub-channels with channel gains given by the diagonal elements of Aj^. 

3) Hybrid D-SVD/P-SVD: A new hybrid D-SVD/P-SVD algorithm is proposed here in the case of multiple 
primary receivers/antennas, i.e., Af,.p > 1. In this algorithm, the secondary user's MIMO channel is first projected 
into the null space of some selected subspace of as opposed to the whole space spanned by in the P-SVD. 
Then, the D-SVD algorithm is applied to this projected channel to obtain transmit precoding matrix as well as 
power allocations that satisfy the interference-power constraints at the primary receivers. There are a couple 
of reasons for using this hybrid D-SVD/P-SVD algorithm. First, it works even if Mt^s < -^r,p under which 
the P-SVD is not implementable. Secondly, in the case of a single-antenna primary receiver, i.e., Mr^p = 1, in 
Theorem [3] and Theorem IH it has be shown that the D-SVD and P-SVD are asymptotically optimal as Pt — > 
and Pt —>■ oo, respectively. It is thus conjectured that when Af,.p > 1, the hybrid algorithm is likely to perform 
superior than both the D-SVD and P-SVD for some moderate values of Pt (as will be verified later by the 
simulation results in Section IVll ). 

One important issue to address for the hybrid D-SVD/P-SVD algorithm is on the selection for the projected 
subspace of G^. Although the optimal selection rule remains unknown, in this paper a heuristic rule is proposed 
as follows. First, each channel gr^ ^ in G is normaUzed by the corresponding ^/Tk so as to make the equivalent 
interference-power constraints at all primary receivers equal to unity. Denote this normalized G as G. Next, let 

1/2 t 1/2 

the SVD of G be expressed as G = QqA ^ [/L,, where the singular values in A^ are ordered by A^, ^ > . . . > 
, = min(Mt M^^p). The hybrid algorithm then projects H into the null space of the space spanned 
by the first b , b < M^, column vectors of U^, denoted by U ^{b), corresponding to the first b largest singular 
values of G, i.e.. 



After applying the precoding matrix V based on the SVD of H^{b), it can be verified that the secondary 
user's MIMO channel is diagonalized and decomposed into min(Mr Mj <j — b) sub-channels. Notice that if 
b < M^, there may be remaining interference at each primary receiver and, hence, the secondary transmitter 
power allocation 5] needs to be adjusted to make the interference-power constraints being satisfied at all primary 
receivers. This can be done by solving a similar problem like Problem P3 in the D-SVD case. 

The rational for the hybrid D-SVD/P-SVD algorithm lies in that only some selected subspace of G^ at the 




(50) 
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primary receivers is removed by the secondary transmit precoding, while the remaining subspace of is 
preserved for the secondary transmission. Notice that the nulled subspace on the average contributes to the most 
amount of interference powers at the primary receivers because it corresponds to the first h largest singular 
values of G. For typical wireless channels, the channel from the secondary transmitter to each primary receive 
anteima (represented by the corresponding row vector of G) is usually subject to variation or fading over time. 
Furthermore, these chaimels may exhibit quite different average chaimel gains (e.g., a near-far situation), and 
some of them may have certain correlations (e.g., for the receive anteimas at the same primary receiver or for 
the primary receives located in the same vicinity). All these factors can result in a more spread distribution 

1/2 

for the singular values in A / , which thereby makes the subspace projection more effective in removing the 
interferences at the primary receivers. Because G is usually independent of the secondary MIMO chaimel H, 
it is unlikely that the nulled subspace of G^ are strongly correlated with H, a phenomenon termed interference 
diversity. Therefore, by selecting a proper value of b, the secondary user has a more flexibility in balancing 
between spacial multiplexing for the secondary transmission and interference avoidance at the primary receivers. 
Notice that the interference-diversity effect is in principle very similar to the multiuser diversity effect previously 
reported in the literature (e.g., [33]-[34]), which is exploited for multiuser rate scheduling in a mobile wireless 
network. 

V. Multi-Channel Transmission 

So far, we have considered the single-channel transmission for both primary and secondary users, and shown that 
even when some primary users are active for transmission, the secondary user is still able to achieve opportunistic 
spectrum sharing with active primary users by utiUzing multiple transmit anteimas and properly designing its 
transmit spatial spectrum. In this section, a more general multi-channel transmission is studied where both 
primary and secondary users transmit over parallel single-channels. This scenario is applicable when, e.g., both 
primary and secondary users transmit over multi-tone or orthogonal-frequency-division-multiplexing (OFDM) 
channels, or alternatively, over consecutive block-fading channels. For such scenarios, in order to achieve optimum 
spectrum sharing, the secondary user needs to first detect the primary user's transmission activities in all available 
dimensions of space, time and frequency and then adapts its transmit resources such as power, rate, and spatial 
spectrum at all of these dimensions. For convenience, in this section we consider the case of one single-antenna 
primary receiver in the CR network although the generalization to multiple primary receivers with single or 
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multiple antennas can also be done similarly like in Section ITVl 

Let Hj, Qj and Sj be defined the same as in Section |lll] but now corresponding to one particular sub-channel 
j,j = 1, . . . ,N, where N is the total number of sub-channels for the multi-band transmission. Each sub-channel 
can be considered as, e.g., one OFDM tone in frequency domain or one fading block in time domain. It is assumed 
that the interference-power constraint 7 is identical at all sub-channels. Problem P2 can be then extended to the 
multi-channel transmission as (P6) 



(51) 



N 

McLximize ^^log2 I + HjSjH^j 

N 

Subject to ^Tr(S'j) < Pt, (52) 

g,Sjg]<jyj, (53) 

SjhO,yj. (54) 

This problem maximizes the total transmit rate over all of N sub-channels for the secondary transmission under a 
total transmit-power constraint and a set of interference-power constraints each for one of N sub-channels. Like 
Problem P2, the above problem can be shown to be convex and, hence, can be solved using standard convex 
optimization techniques. In the following, the Lagrange dual-decomposition method is applied to decompose 
Problem P6 into N subproblems all having an identical structure. The main advantage of this method lies in that 
if the same computational routine can be simultaneously applied to all N subproblems, the overall computational 
time is maintained regardless of A^. Other applications of the Lagrange dual-decomposition method for resource 
allocation in communication systems can be found in, e.g., [35]-[43]. 

For the problem at hand, the first step is to introduce the non-negative Lagrange multiplier u associated with 



the power constraint in (I52b and to write the Lagrangian of the original (primal) problem as 

N / ^ \ 

C{{S,},u) = J^log2 \i + H,S,h]\ - y \y^Tr{Sj)-Pt . (55) 

Let Vj be the set specified by the remaining constraints in ([53]) and ([54]) corresponding to sub-channel j, j = 
1, . . . ,N. The Lagrange dual function is then defined as 

/(z.)= max C{{Sj},u). (56) 

The dual problem is then defined as min,y>o f{i^)- It can be verified that the strong duality holds for the problem 
at hand, and hence the duality gap is zero. Therefore, the primal problem can be solved by first maximizing the 
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Lagrangian £ to obtain the dual function /(i^), and then minimizing f{ij) over all non-negative values of u. 

Consider first the problem for obtaining f{v) with some given u. It is interesting to observe that /{u) can be 
also written as 

N 

where 

f'{u) = max log2 I + HjSjH] - uTr{Sj), j = 1, . . . , iV. (58) 



Hence, f{v) can be obtained by solving N independent subproblems, each for f'j{i^), j = l,...,N. Each 
subproblem has the same structure and hence can be solved by using the same computational routine. This 
practice is usually referred to as the dual decomposition. From (1581 ). each of these subproblems can be written 
more explicitly as (P7) 

rt 



Maximize log2 



I + HjS,H. 



i/Tr(5j) (59) 



Subject to 9jSjg] < 7, (60) 
Sj h 0. (61) 

Compared with Problem P2, the above problem has the difference in that the transmit-power constraint in ([TTI) 
is removed and a new term z^Tr(S'j) is subtracted in the objective function. This is a consequence of a total- 
power constraint over all N sub-channels in the multi-channel case instead of the per-sub-channel-based power 
constraint in the single-channel case. It can be verified that Problem P7 is also convex and hence can be solved 
by using, e.g., the SOCP or SDP if the secondary user's channel is MISO (beamforming is optimal), or the 
interior-point method if the secondary user's channel is MIMO (spatial multiplexing is optimal). 

The dual variable u can be considered as a common price applied to all sub-channels for regulating their 
allocated transmit powers. It can be verified that in Problem P7, for each sub-channel, a larger u will result in 
a smaller power consumption Tr(S'j) and vice versa. Hence, the optimal u can be found by a bisection search 
by comparing the optimal sum-power over all N sub-channels for a given v with the transmit-power constraint 
Pt. In summary, the following algorithm (A2) can be used to solve Problem P6: 

• Given G [0, ;>] 

• Initialize fmin = 0, fmax = 
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• Repeat 

1. Set V ^ ^(t'mm + t'max)- 

2. Solve Problem P7 for j = 1, . . . , independently to obtain {Sj}. 

3. If Y.j=\ < Pt, set fmax ^ otherwise, 

• Until J^max — ^xmn < (J;/ where bp is a small positive constant that controls the algorithm accuracy. 
Algorithm A2 can also be applied when the D-SVD or P-SVD algorithm is used to approximately solve each 

subproblem P7. The selection of the D-SVD or P-SVD at each sub-channel depends on the dual variable v and 
the channels Hj and g^. Because Hj and may change over different sub-channels, in order to exploit this 
diversity, the best selected SVD-based algorithm (the D-SVD or P-SVD) may also change from one sub-channel 
to another. 

VI. Simulation Results 

For the simulations, it is assumed that in the CR network both the channels from the secondary transmitter to 
the primary receivers G, and to the secondary receiver H, are drawn from the family of random vectors/matrices 
for which each element is generated by independent CSCG variable distributed as CAA(0, 0.1) for G, and CAA(0, 1) 
for H, respectively. The secondary user's transmit-power constraint Pt sweeps from 1 to 100, which is equivalent 
to a range of average SNRs per receive antenna (measured when the actual transmit power P is equal to Pt) from 
dB to 20 dB. Monte Carlo simulations with 1000 randomly generated channel-pairs (G, H) are implemented, 
and the average achievable rates are plotted versus the SNR values. The results are shown for the following cases: 

A. Capacity With (w/) versus Without (w/o) Interference-Power Constraint 

Fig j2] compares the capacity for the secondary user's transmission w/ and w/o the interference-power constraint 
at a single-antenna primary receiver. In the case that the interference-power constraint is applied, 7 is equal to 
0.01. For the secondary receiver, it is assumed that M^^s = 1> and two antenna configurations are considered for 
the secondary transmitter: (1) Single-input single-output (SISO) case: Mt^s = 1; (2) MISO case: Mt^s = 4. For 
both SISO and MISO secondary user's channels, it is observed that the interference-power constraint reduces the 
capacity of the secondary transmission, especially at the high-SNR regime. However, the capacity improvement 
by adding more transmit antennas at the secondary transmitter is observed to be substantial. This is so because in 
the SISO case, as SNR increases, the capacity is eventually limited by the interference -power constraint instead of 
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the secondary user's own transmit power constraint; while in the MISO case, the secondary transmitter is able to 
apply transmit beamforming so as to avoid the interference at the primary receiver and at the same time, guarantee 
its own rate increase with SNR. Therefore, in addition to antenna array and diversity gains in the conventional 
MISO point-to-point transmission, multi-antennas at the secondary transmitter play another essential role in a 
CR network, i.e., adapting transmit spatial spectrum away from the interfering direction to the primary receiver, 
so as to achieve a high spectral efficiency for the secondary transmission. 

B. D-SVD versus P-SVD 

Fig. [3] shows the achievable rates of two SVD-based algorithms, the D-SVD and P-SVD, and compares them 
to the capacity w/ and w/o the interference-power constraint, as well as the achievable rate by a simple design 
for the secondary transmit spatial spectrum, termed white spectrum, for which equal power is allocated to all 
transmit antennas, i.e., S = jj^I- It is assumed that there is a single-antenna primary receiver, 7 = 0.1, and 
Mt^s = Mr^s = 2. It is observed that the achievable rates by the D-SVD and P-SVD are close to the capacity 
w/ interference -power constraint at the low- and high-SNR regime, respectively. At the high-SNR regime, both 
the capacity by the optimal S obtained by numerical algorithm and the achievable rate by the P-SVD increase 
linearly with log2 SNR by a factor of m.m{Mr^s, Mt^s — 1) = min(2,2 — 1) = 1 while the achievable rate by 
the D-SVD eventually gets saturated, both in accordance with Theorem |4l Substantial rate improvements by the 
optimal as well as the D-SVD/P-SVD solutions over the simple white spectrum are also observed from low to 
high SNR values. 

C. Hybrid D-SVD/P-SVD 

Fig. |4] compares the achievable rates of the hybrid D-SVD/P-SVD with different values of b. The capacity 
w/ and w/o the interference-power constraint are also shown for comparison. It is assumed that the number of 
primary receivers K = 2 each with a single receive antenna, F^, = 0.1, A; = 1, 2, and Mt^s = -^r,s = 4. Because 
there are two single-antenna primary receivers, the hybrid algorithm can choose to project the secondary user's 
channel H into either the whole space of by taking 6 = 2 (or the P-SVD) or some selected subspace of 
with 5 = 1, or simply choose not to project at all with 6 = (or the D-SVD). It is observed that as SNR increases, 
the optimal SVD-based algorithm changes from the D-SVD to the hybrid D-SVD/P-SVD with 6 = 1, and finally 
to the P-SVD. This result is consistent with the conjecture we have made previously in Section IIV-B.3I The 
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rational here is that as the transmit power increases, the secondary user's transmit spatial spectrum also needs to 
allocate a larger portion of total transmit power into the projected channel in order not to violate the prescribed 
interference-power constraint at the primary receiver. 

Fig. [5] compares the achievable rates of the hybrid D-SVD/P-SVD with the D-SVD and P-SVD as the number 
of single-antenna primary receivers K increases from 2 to 10. The interference-power constraint = 0.01 is 
assumed to be equal at all primary receivers, and Pt is fixed as 10, and ^ = M^^s = 4. Notice that for each 
given K, b for the hybrid algorithm can take values from to min(Mt ,5, -?^)- If 6 = 0, the hybrid algorithm 
becomes the D-SVD; while if 6 = min(Mt K), it becomes the P-SVD. In this figure, the achievable rate shown 
for the hybrid algorithm for each K is obtained by taking the maximum rate over all possible values of b. Notice 
that if > Mt,s, the achievable rate by the P-SVD becomes zero because the ZF-based channel projection 
results in a null projected channel H ±^ in ( |49l ). It is observed that as K increases, the achievable rate by any 
of these three proposed algorithms decreases due to more added interference-power constraints. However, the 
hybrid algorithm by dynamically selecting b achieves substantial rate improvements compared with the D-SVD 
and P-SVD because it effectively exploits the interference-diversity gain in the CR network. As K becomes much 
larger than Mt^s, it is observed that the achievable rates by the hybrid algorithm and the D-SVD both converge 
to be identical. However, such convergence is observed to be very slow with respect to K. 

D. Multi-Channel Transmission 

A multi-tone or OFDM-based broadband transmission system is considered for a secondary user with Mt^s = 
Mr,s = 2 and a single-antenna primary receiver. The OFDM channels for both primary and secondary trans- 
missions are assumed to have = 64 tones, and four equal-energy, independent, and consecutive multi-path 
delays. The interference-power constraint at the primary receiver is set to be 7 = 0.1 at all tones. In this case, 
the secondary user is able to allocate variable transmit rate, power, and spatial spectrum at different OFDM tones 
based on the channels Hj and g^, j = I, . . . ,N. Fig. [6] shows the secondary user's capacity w/ and w/o the 
interference-power constraint, and the achievable rate of the per-tone-based suboptimal algorithm that selects the 
best SVD-based algorithm (the D-SVD or P-SVD) at each tone. It is observed that by exploiting the frequency- 
selective fading, the per-tone-based suboptimal algorithm achieves smaller gaps from the actual capacity w/ 
interference-power constraints compared with the case of single-channel transmission previously shown in Figl3] 
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VII. Concluding Remarks and Future Directions 

Transmit optimization for a single secondary MIMO/MISO link in a CR network under constraint of oppor- 
tunistic spectrum sharing is considered. The capacity of the secondary link is studied under both the secondary 
transmit-power constraint and a set of interference-power constraints at multiple primary receivers. Multi-antennas 
are exploited at the secondary transmitter to optimally tradeoff between throughput maximization and interference 
avoidance. In the case of MISO secondary user's channel, beamforming is shown to be the optimal strategy for 
the secondary transmitter. In the case of MIMO secondary user's channel, the capacity-achieving transmit spatial 
spectrum under the interference-power constraint in general does not diagonalize the MIMO channel and hence 
requires sophisticated encoding and decoding methods. Two suboptimal algorithms based on the SVD of the 
secondary user's MIMO channel, namely the D-SVD and the P-SVD, are proposed to tradeoff between capacity 
and complexity. In the case of multiple primary receivers/antennas, a hybrid D-SVD/P-SVD algorithm is also 
proposed to exploit the inherent interference diversity gain in a CR network. The developed algorithms are also 
extended to the multi-channel transmission whereby the secondary user is able to employ transmit adaptations 
in all space, time and frequency domains so as to achieve optimum opportunistic spectrum sharing. 

An interesting extension of this paper is to consider more generahzed interference constraints at each primary 
receiver, especially when it is equipped with multiple receive antennas. Instead of considering either the per- 
antenna-based power constraint or the total-power constraint over all receive antennas, the interference constraint 
at each primary receiver can also be more specifically related to the primary user's channel capacity. More 
interestingly, both primary and secondary transmit adaptations can be jointly optimized under mutual interferences 
given their unequal priorities for spectrum utilization, through either a centralized control or distributed self- 
adaptations based on the principle of game theory. Results and algorithms developed in this paper can also be 
extended to the case where only statistical knowledge on the channel, instead of instantaneous channel knowledge 
as assumed in this paper, is available at the secondary transmitter for optimum design of transmit spatial spectrum 
under some long-term average interference-power constraints. Furthermore, although this paper addresses a single 
pair of secondary transmitter and receiver, it can also be extended to incorporate multiple secondary users that 
jointly share transmit spectrum with the primary users by considering different models of the secondary network, 
e.g., the uplink MIMO multiple-access channel (MAC), the downUnk MIMO broadcast channel (BC), or the 
distributed MIMO interference channel (IC). 
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Appendix I 
Proof of Theorem [T] 

Consider the k-th primary transmission represented by 

yi,{n) = HkXk{n) + Zk{n) + qi,{n), (62) 

where G £.MkxNk denotes the k-th primary user's channel, and are the number of antennas at its 
receiver and transmitter, respectively; 2/fc(n) and Xk{n) are the received and transmitted signal vector, respectively; 
qj^{n) is the additive noise vector at the receiver that assumed to be distributed as Zk{n) CJ\f{0, (pkl)', Qfc('^) 
is the interference from the secondary transmitter expressed as 



Qkin) = Gkx{n). 



(63) 



It is easy to verify that the covariance matrix for qk{n), denoted by Qj^ = E[g^(n)g^(n)], is equal to GkSG\, 
where S is the transmit covariance matrix for the secondary user. Let the transmit covariance matrix of the 
k-ih primary user be denoted by S^, = E[a;fc(n)a;|.(n)]. The capacity of the k-ih primary transmission (in 
bits/complex dimension) can be then expressed as 

Ci = log2 |/ + {(Pkl + Qfc)"'/' HkSkHl {(t^kl + Qkr^'^\ ■ (64) 

On the other hand, if the secondary transmitter if off, i.e., the interference from the secondary transmitter ^^(ra) 

is nonexistent, the capacity of the fc-th primary transmission is expressed as 

1 



Co = log 



0fc 



(65) 



The capacity loss of the A;-th primary transmission due to the secondary transmission is then equal to C2 — Ci. 
In order to find an upper-bound for such capacity loss, the following equalities/ineuqalities are provided for Ci. 



Ci 



log2 I + S]/^hI{<PuI + Qkr^ H^S]!^ 



> log2 

= log2 

> log2 



/ + 
/ + 



1 



1 



(66) 
(67) 

(68) 
(69) 
(70) 

where (1661 ) and (1681 ) are due to log |/ + AB\ = log |/ + BA\, (167] ) is from the facts that Q/. ^ r^I and 
log |/ + AXiA^l > log \I + AX2A'f\, if ^1^X2^ 0. Using the proof is completed. 



min(Mfc,Arfc)log2 ( 1 + ^ ) + C2, 
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Appendix II 
Proof of Lemma [T] 

Because Problem P2 is convex, the optimal S must satisfy the Karush-Kuhn-Tacker (KKT) conditions [23] 
written at below: 

h^(l + hSh^^ \ + * = iig^g + vl, (71) 
v{Tr{S)-Pt) = 0, (72) 

^Ji{gSg^--^) = o, (73) 

Tr($S) = 0, (74) 

where v, ^, and $ are the Lagrange multipliers associated with the constraint ([TT]) . ([T2l) and (fT3l ). respectively, 
I/, /X > 0, $ ^ 0. First, consider the case of u = 0, from (1721 ). it follows that P < P^, i.e., P is not limited by the 
secondary transmit-power constraint, but instead, by the interference-power constraint. From dTT] ) and because 
$ ^ 0, it follows that g must be in parallel to h, i.e., g = ch, where c is a constant. In this case, it can be 
easily shown that the optimal S is S = j^^^^j-^h'^h, i.e., Rank(S) = 1. Next, consider the case of v > 0, i.e., 
P = Pf. In this case, the RHS of (TtT] ) has a full rank of ^ regardless of fi. It then follows that at the LHS of 
(TtT]) . since the first term has a unit rank, Rank($) > Mt^s — 1- Since S ^ and $ ^ 0, from (1741 1 it follows 
that Rank(S) + Rank($) < Mt^s- Hence, Rank(S) < 1, and the proof is completed. 

Appendix III 
Proof of Lemma [2] 

Let the optimal beamforming vector v' = a^'g + (3v'b, where 1) g = 0. It can be then verified that replacing h 
with /ix does not increase the interference power Hfl't'lP in (IT6l ) since {h±)'^g = 0, but always helps to improve 
the SNR at the secondary receiver in ([T4b since > The proof is then completed by the 

fact that the secondary user's capacity always increases with the receiver SNR. 
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Appendix IV 
Proof for Optimality of P-SVD when 7 = 



We show that if in the interference-power constraint (112] ) of Problem P2, 7 = 0, the P-SVD is indeed optimal. 
First, using the constraint gSg^ = 7 = 0, ([TOl l can be simplified as 



log2 
log2 
log2 



I + HSH^ 
I + {H^ + Hgg^)S{H^ + Hgg^)^ 
I + H^S{H^)^ 



(75) 
(76) 
(77) 



where (1761 ) is from (|27] ). Using (TTTl) . Problem P2 with 7 = can be rewritten as 



Maximize 
Subject to 



log2 I + H^S{H^)^ 
Tr(5) < Pu 
gSg^ = 0, 
5^0. 



(78) 
(79) 
(80) 
(81) 



Notice that without the constraint (1801 ). the optimal S* for the above problem should have the same form of 
U±'S(U±)'^ as by the P-SVD. At last, it remains to check whether gS*g^ = 0. This is so because g'^U ±_ = 0. 

Appendix V 
Proof of Lemma [3] 



The following equalities/inequalities hold: 



log2 



I + HSH 



log2 I + {H^ + Hgg^)S{H^ + Hgg 



< log2 

< log2 



/ + 2H^S{H^)^ + 2Hgg^S{Hgg)^ 



I + 2H^S{H_ 



+ log2 



I + 2Hgg^S{Hgg) 



(82) 
(83) 
(84) 
(85) 



where (ES is from (EH), dH is by {A + B)S{A + B)t ^ 2{ASA^ + BSB^) for 5 ^ 0, ([85]) is by log \I + 
51 + 521 < log |/ + Sil + log \I + S2\ for 5i, S2 ^ 0. Let Vs denote the set for S specified by ([12]) and 
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fT3] ). it then follows that 



Ropt = max log2 

S&Vs 



I + HSH^ 



< max log2 

S&Vs 

< max log2 
SeVs 



(86) 

I + 2H^S{H^)^\ + 0{1), (87) 
2(^I + H^S{H^)^^\ + 0{1), (88) 
< i?P-svD + Af,,, + 0(l), (89) 



where (l87l) follows from ([85] ) by noticing that due to (fT2l ) the second term on the RHS of (1851 ) is upper-bounded by 
some finite constant. From ( 1891 ) and because i?opt > ^p-svd> it follows that limp^^^oo iH^Tf ~ limPt^oo log ■ 
Since it has already been verified in the main text that limp^^oo ^og^K ~ ^^^i^r,s, Mt^s — 1)> the proof is 
completed. 
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Fig. 1. A cognitive radio (CR) network where the secondary user shares the same transmit spectrum with K primary users. 
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Capacity comparison for the secondary transmission w/ and w/o the interference-power constraint at a single-antenna primary 
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Fig. 3. Comparison of the achievable rates for the secondary transmission by the D-SVD and P-SVD in the case of a single-antenna 
primary receiver, and Mt,s = Mr,s = 2. 
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Fig. 4. Comparison of the achievable rates for the secondary transmission by the hybrid D-SVD/P-SVD with different values of b, in 
the case of two single-antenna primary receivers, and Mt,s = Mr,s = 4. 
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Fig. 5. Comparison of the acliievable rates for the secondary transmission by the hybrid D-SVD/P-SVD, the D-SVD and the P-SVD 
for different number of single-antenna primary receivers K, in the case of Mt,s = Mr,s = 4 and a constant secondary power constraint 
Pt = 10. 
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Fig. 6. Comparison of the achievable rates for the secondary transmission in a multi-tone-based broadband system with iV = 64, a 
single-antenna primary receiver, and Mt,s = Mr,s = 2. 



